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Abstract 

In this paper we study the optimal placement and optimal number of active relay 
nodes through the traffic density in mobile sensor ad-hoc networks. We consider a 
setting in which a set of mobile sensor sources is creating data and a set of mobile 
sensor destinations receiving that data. We make the assumption that the network is 
massively dense, i.e., there are so many sources, destinations, and relay nodes, that 
it is best to describe the network in terms of macroscopic parameters, such as their 
spatial density, rather than in terms of microscopic parameters, such as their individual 
placements. 

We focus on a particular physical layer model that is characterized by the follow- 
ing assumptions: i) the nodes must only transport the data from the sources to the 
destinations, and do not need to sense the data at the sources, or deliver them at the 
destinations once the data arrive at their physical locations, and ii) the nodes have 
limited bandwidth available to them, but they use it optimally to locally achieve the 
network capacity. 

In this setting, the optimal distribution of nodes induces a traffic density that 
resembles the electric displacement that will be created if we substitute the sources and 
destinations with positive and negative charges respectively. The analogy between the 
two settings is very tight and have a direct interpretation in wireless sensor networks. 



1 Introduction 

Various approaches inspired by physics have been proposed to deal with the routing problem 
in massively dense wireless sensor networks. Starting with the pioneering work of Jacquet 
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(see [H [2]) who used ideas from geometrical optics to deal with the case of one source and 
one destination and a distribution of relay nodes. 

Approaches based on electrostatics have been studied in [31 HI El E] (see the survey [3| and 
references therein) to deal with the case of a distribution of sources and destinations with a 
density of relay nodes in a static environment. 

The development of the theory of routing in massively dense wireless sensor networks has 
emerged in a complete independent way of the theory developed within the community of 
road traffic engineers, introduced in 1952 by Wardrop [11] and Beckmann [12], see also [131 
pp. 644, footnote 3] for the abundant literature of the early 50's, and which is still an active 
research area among that community, see [HI UHl HH [TTl [18] and references therein. 

Only very recently, new approaches based on road traffic theory have been studied in [81191 [TO] 
defining the Wardrop equilibrium and a characterization of the Wardrop equilibrium in this 
type of setting and including a geometrical characterization of the flow of information for 
some particular cost functions. 

Consider that a spatially distributed set of mobile sources is creating data that must be 
delivered to a spatially distributed set of mobile destinations. In this context, our objective 
is to study the optimal traffic distribution and to find the minimum amount of relay nodes 
needed to transport the data from the sources to the destinations. 

The main contribution of this work is to address this problem for a mobile context where 
we analyze the cases where (i) only sources are mobile and the destinations are static as it 
would be in the case when the aggregation centers are fixed and the sensor nodes have the 
capability to move, (ii) the case when both sources and destinations are mobile, and given 
that the mathematics involve are similar (iii) we also analyze the case when the sources are 
static and the destinations are mobiles. 

In section 12.11 we first analyze the case where the sources and destinations are static and 
spatially distributed in the one dimensional line in order to illustrate the behavior of the 
fiow function and the optimal location of the relay nodes on this simple case. We analyze 
in section 12.21 the case when the sources and destinations are static and spatially distributed 
in the two dimensional plane. In section [3] we analyze the case where the sources and 
destinations can move with a deterministic velocity and we are able to find the optimal fiow 
of information and to give the optimal spatial density of the relay nodes at each period of 
time. Within this section we also give an example with numerical results related to the 
previous mobility setting. In section HI we are able to find the optimal spatial density of the 
relay nodes for the Brownian mobility model. In section [6] we summarize the main results 
and future perspectives for the continuation of our work. 
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2 The Model 



We first consider the one dimensional case in order to explain the main concepts involved in 
our model and how this concepts can be extended to the two dimensional case in order to 
obtain the optimal deployment of the relay nodes in a wireless sensor network. 

2.1 Fluid Approximations 

Consider a grid area network that contains wireless sensor nodes. As a first approach we 
consider the line segment [0, L], which will be the geographical reference of a network. 

We consider the continuous node density function ri{x), measured in nodes/m, such that the 
total number of nodes on a segment [^0)^i]) denoted by N{io,ii), is 



We consider as well the continuous information density function p{x) , measured in bit /(s ■ m) , 
generated by the sensor nodes such that 

• At locations x where p{x) > there is a fraction of data created by the sensor sources, 
such that the rate with which information is created in an infinitesimal area of size de, 
centered at position x, is equal to p(x) de. 

• Similarly, at locations x where p{x) < there is a fraction of data received at the sensor 
destinations such that the rate with which information is received by an infinitesimal 
area of size de, centered at position x, is equal to —p{x) de. 

We assume that the total rate at which sensor destinations have to receive data is the same 
as the total rate which the data is created at the sensor sources, i.e., 



Notice that if we have an estimation of the proportion of packet loss in the network, we can 
ponderate our function p in order to adequate it to equation ([1]). 

Consider the continuously differentiable traffic flow function T{x), measured in bps/m, such 
that its direction (positive or negative) coincides with the direction of the fiow of information 
at point X and |T(x)| is the rate at which information propagates at position x, i.e., \T{x)\ 
gives the total amount of traffic that is passing through the position x. 




L 
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Next we present the flow conservation condition. For information to be conserved over a 
segment [£o,^i], it is necessary tliat tlie rate witli wliicli information is created over the 
segment, is equal to the rate with which information is leaving the line segment, i.e., 



T(£i)-T(4) = / p{x)dx 

Jin 



The integral on the right hand side is equal to the quantity of information generated (if 
it's positive) or demanded (if it's negative) by the fraction of sensor nodes over the line 
segment [^c^i]- The expression T(£i) — T(£o), measured in bps/m, is equal at the rate with 
which information is leaving (if it's positive) or entering (if it's negative) the segment [io, ii]. 
This holding for any line segment, it follows that necessarily. 

The problem considered is to minimize the number of nodes A^(0, L) in the line segment [0, L], 
needed to support the information created by the sensor sources and received by the sensor 
destinations subject to the flow conservation condition given by equation ([2]) and imposing 
that there is no flow of information leaving the network, i.e., T{0) = and T{L) = 0. Thus 
the system of equations that model our problem in the one-dimensional case is given by the 
system of equations: 

MmN{0,L) = / ri{x)dx, (3) 



subject to — ; — = p{x) in(0,L), (4) 

dx 

T(0) = and T(L) = 0. (5) 

Notice that in the one- dimensional case, there is no minimization problem because only by 
using the constraints @ and ([5]), we obtain only one solution. As we will see, this will not 
be the case for the two-dimensional case. 

We suppose that the proportion of sensor nodes r]{x) in an area of infinitesimal size de, 
centered at location x, needed as relay nodes, will be proportional to the traffic flow of 
information that is passing through that region, i.e., ri{x)de = \T{x)\°'de where a > is 
a fixed number called the relay-trafEc constant. Then the optimal placement of the relay 
nodes in the network will be given by 7]*{x) = |T*(x)|", where the traffic flow function T*{x) 
is the optimal traffic flow function, given by the solution of the previous system of equations. 
Furthermore, the optimal total number of relay nodes A^*(0, L) needed to support the optimal 
traffic flow function T*{x) in the network will be N*{io,ii) = J^^ r]{x) dx = /;^'|T(x)|" rfx. 
Let us see an example to illustrate the previous framework. 

Example 1.- Suppose that we can divide the line segment [0, L] in two parts: (i) in the flrst 
part [0,L/2] there will be a uniform information density function generated by the sensor 
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sources, given by p{x) = 1 bps/m and (ii) in the second half [L/2, L] there will be a uniform 
information density function received at the sensor destinations given by p{x) = —1 bps/m 
(See Figure 1). From the equations (jlj) and ([5]) we obtain that the optimal traffic flow 
function will be given by T*{x) = x bps/m for all x G [0,L/2] and T*{x) = L — x bps/m 
for all X G [L/2,L] with positive direction (See Figure 2). If we assume that the relay- 
traffic constant a = 2, then the optimal placement of the relay nodes needed to relay the 
information from the sources to the destinations on the network will be given by ri*{x) = 
for all X G [0, L/2] and T*{x) = (L - xf bps/m for all x G [L/2, L]. 

The optimal total number of relay nodes A^*(0, L) needed to support the optimal traffic T*{x) 
will be given by N{L) = J^^^ x^ dx + j^JL - xf dx = L^/12. 
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2.2 The two dimensional case 

Consider a grid area network D in the two dimensional plan^ X x Y the continuous infor- 
mation density function p(x), measured in bps/m^, such that at locations x where p(x) > 0, 
there is a distributed data created by sensor sources, such that the rate with which infor- 
mation is created in an infinitesimal area of size dA^, centered at location x, is p(x) dA^. 
Similarly, at locations x where p(x) < 0, there is a distributed data received at sensor des- 
tinations, such that the rate with which information can be treated by an infinitesimal area 
of size dA,,, centered at location x, is equal to — p(x) dAf,. 

The total rate at which sensor destinations must process data is the same as the total rate 
which the data is created at the sensor sources, i.e., 

J p(x) dx = 0. 

XxY 

Consider the continuous node density function r7(x), measured in nodes/m^, defined so that 
the number of relay nodes in an area of infinitesimal size dA^, centered at x, is equal to 
rjix.) dA^. 

The total number of nodes on a region A, denoted by N{A), is then given by 

N{A) = j ri{x) dx. 

A 

Consider the continuous traffic flow function T(x), measured in bps/m, such that its direction 
coincides with the direction of the fiow of information at point x, ancl§ ||T(x)|| is the rate 

^We will denote with bold fonts the vectors and x = (x, y) will denote a location in the two dimensional 
space XxY. 

^ The norm ||-|| is the Euclidean norm, i.e., for a vector x — {x,y), its norm will be ||x|| = y'x^ + y^. 
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with which information rate crosses a hnear segment perpendicular to T(x) centered on 
X, i.e., ||T(x)|| e gives the total amount of traffic crossing a linear segment of infinitesimal 
length e, centered at location x, and placed vertically to T(x). 

Next we present the fiow conservation condition (see e.g. [HIH], for more details about this 
type of condition). For information to be conserved over a domain D of arbitrary shape 
on the X X Y plane, with smooth boundary dD, it is necessary that the rate with which 
information is created in the area is equal to the rate with which information is leaving the 
area, i.e., 

J p(x) dD = ^[T ■n{x)]di 

D dD 

The integral on the left is the surface integral of p(x) over the domain D. The integral on 
the right is the path integral of the inner product T ■ n over the boundary dD. The vector 
n(x) is the unit normal vector to dD at the boundary point x G dD and pointing outwards. 
Then the function T-n (x), measured in bps/m^, is equal at the rate with which information 
is leaving the domain D at the boundary point x. 

This holding for any (smooth) domain D, it follows that necessarily 

_ , 9Ti(x) 9T2(x) 

where "V-" is the divergence operator. 

Thus the problem considered is to minimize the quantity of nodes N{D) in the grid area 
network D needed to support the information created by the distribution of sources subject 
to the fiow conservation condition, i.e., our problem is given by the system of equations: 

Mm N{D) (7) 
subject to V ■ T = p(x). (8) 

Toumpis and Tassiulas in [3] focus on a particular physical layer model characterized by the 
following assumption: 

Assumption : A location x where the node density is //(x) can support any traffic fiow vector 
with a magnitude less or equal to a bound ||T(x)||inax which is proportional to the square 
root of the density, i.e. ||T(x)|| < ||T(x)||ma2; = K^J d{x.). 

The validity of Assumption 1 depends on the physical layer and the medium access control 
protocol used by the network. Although it is not generally true, it holds in many different 
settings of interest. For example, in [3J Toumpis and Tassiulas give an example of network 
where nodes are placed in a perfect square grid of m x m nodes and each node can listen 
to transmissions from its four nearest neighbors. They give a simple time division routine so 
that the network of nodes can support a traffic on the order of m. As another example, 
in [20] it was shown that the traffic that can be supported in the above network, if nodes 
access the channel by use of slotted Aloha instead of time division, is Tiocai = K x W x m. 
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where nodes transmit data with a fixed global rate of W bps, is a constant-smaller than 
1/3-that captures the efficiency of Aloha. Finally, in [19] it was shown that a network of n 
randomly placed nodes can support an aggregate traffic on the order of ^/n/~\ogn under 
a more realistic interference model that accounts for interference coming from arbitrarily 
distant nodes. The logarithm in the denominator appears due to the proving methodology 
of [19] , and it has been shown [21] that it can be dispensed off, by use of percolation theory 
in the proofs. 

Tassiulas and Toumpis prove in [3] that among all traffic ffow functions that satisfy V • T = 
p, the one that minimizes the number of nodes needed to support the network, must be 
irrotational, i.e., 

V X T = 0. (9) 

where "Vx" is the curl operator. 



3 Model in motion 

In our work we do a parallel to Electromagnetism by considering moving sensor nodes or 
moving distribution of sensor nodes. In that sense we ffrst deffne some terms directly related 
to Electromagnetism. In our model, the information density function p, the traffic ffow 
function T, and the node density function rj, defined previously may depend on time, i.e., 
p = p(x, t), T = T(x, t), and rj = //(x, t).We will consider our problem within a window of 
time t G [ti, tj] where ti is the initial time and tf is the final time. 

We define the continuous node current J as the density of sensor nodes p(x, t) in the position 
X multiplied by the nodes average drift velocity v(x, t), i.e., 

J = p(x,t) v(x,t). 

The rate at which nodes leaves an area (or volume) V, bounded by a curve (or surface) 
S = dV , will be given by 

J • dS (10) 



s 

Since the information density function is conserved in the plane this integral must be equal 
to 

lj.,S^4l,.n,S^-lpV. (11) 

From the divergence theorem and imposing the equality between the equations flTU]) and 
flTT]) . we obtain the equivalent to Kirchhoff's current law: 

Notice as well that 

V ■ J = V ■ (p v) = V ■ Vp + pV ■ V. 
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We assume that we know the initial distribution of the sensor sources and the sensor desti- 
nations at time denoted by po- Thus we obtain the following system of equations: 

(TV\! m+^-'^P + P'^-^ = ^ inDx(0,T) 
^ ^1 p(0)=Po onDx{0}. 

The previous system of equations is known in the partial differential equations literature 
as the linear transport equation with initial condition for which there exists a solution (see 
Proposition 11.1 of [22j). 

Notice that given the initial distribution of the sources and destinations and the velocity of 
the distribution, we are able to compute the evolution of the distribution of sensor sources and 
sensor destinations on time t G [0, T). The velocity of sensor sources and sensor destinations 
may be estimated by having some previous knowledge on the behavior of these sources and 
destinations in our network. 

We want to minimize the number of relay nodes N{D) in the grid area network D needed 
to support the information created by the distribution of sensor sources and received by the 
distribution of sensor destinations subject to the flow conservation condition, and knowing 
that the distribution of mobile sensor nodes and mobile sensor destinations is the solution 
to the system of equations (TE). Thus our problem reads for all t E [0,T] 

Min N{D,t)= / r/(x,t)rfx= = / |T(x,t)|2rfx (12) 
Jd J d 



subject to V ■ T(x, t) = p{'s.,t) in D, (13a) 
T ■ n = on aD. (13b) 
where p(x, t) is the solution to the problem (TE). 

We recall that Tassiulas and Toumpis in [3J proved that among all traffic flow functions that 
satisfy equation fll3al) . the one that minimize the number of nodes needed to support the 
network, must satisfy 

V X T = 0. (14) 

Using Helmholtz's theorem (also known as fundamental theorem of vector calculus) to last 
equation 0141) we obtain that there exists a scalar potential function Lp such that 

-V(/? = T. (15) 

Replacing this function into the conservation equation (I13ap we obtain that 

- = p (16) 

and this holds for all t G [0, T). 
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We impose that no information is leaving the considered domain in equation (]13bp and 
from equation (fT5l) this condition translates into Vv? ■ n = 

From equation ( fT6i) and last condition we obtain the following system 

(LE)( -^^ = ^ ^^^^ (17) 
^ ' y Vv9 • n = on oD. ^ ' 

which is the Laplace equation with Neumann boundary conditions. 

If the function / is square integrable then the Laplace equation with Neumann boundary 
conditions has a unique solution in H^{D)/M.. 

In summary in order to solve our problem given by the equations ( IT2|) . fll3al) . (]13b[) . and 
(TE) we need to 

1. Solve the system of equations (TE), 

2. Put the solution as input into the system of equations (LE), 

3. Solve the system of equations (LE). 

Even as it looks complicated system we give an example where you can get explicit solutions. 

Example 2.- Due to presentation effects we consider the one dimensional case during T = 2 
hours. 

We consider an initial distribution of sensor sources Pq and an initial distribution of sensor 
destinations pg on the positive real line [0, +oo), and we scale them to be probability distri- 
butions so it represents in each location the proportion of sensor sources or the proportion 
of sensor destinations respectively: 

p+ = k,e-^-^^" and p- = -k^e-^^-''^". 

where ki and k2 are normalization factors given by ki = 0Perfc(-3) ' ^2 = 0fcrfc(-io) ' ^here 
erfc(x) is the complementary error function defined as erfc(x) = J^°° e~'^^ ds. 

We consider that the nodes average drift velocity is given by v{x,t) = x. We can think of 
a highway where the cars are equipped with sensors and while they are advancing on the 
highway they can go faster and faster. 

[1]- We first need to solve the transportation equation system (TE), that in our example 
reads 

t + ^ = onM+x[0,T) 
p(0) = p++Po onM+ 

Using the method of characteristics we obtain that the information density function over 
time is given by p{x,t) = p'^{x,t) + p~{x,t) where 

p+{x,t) = A;ie-(^^"'-3)'-* 
10 



Figure 1: Distribution of the sources and destinations in the same hne 



The solution combining the sources and destinations information density function over time 
are showed in Fig. [TJ 

O- We put the solution as input into the system of equations: from the conservation equation 
we obtain 

dT{x,t) _ dT+{x,t) dT-{x,t) 
dx dx dx 

dT+{x,t) + dT-{x,t) 

where = p and = — p . 

ox ox 

with initial condition that the flow is zero at the boundary point zero, i.e. T{0,t) = 0. 

O- We solve the laplacean system of equations: Then the optimal traffic flow function is 
given by T*{x, t) = T~^{x, t) + T~{x, t) where 

T-{x,t) = -j;k2e-^^'''-''^'-'dx. 

Thus the minimal number of active relay nodes needed to support the optimal flow at every 
time t will be given by 

r+oo 

N*{t) = / \T*{x,t)fdx 
Jo 

which can be solved numerically. 

In next section we present another type of mobility model where we consider the randomness 
in the mobility of the users. 



4 Brownian Mobility Model 

One of the most used mobility models used in networks is the Random Walk Mobility Model 
also known as the Brownian Mobility Model (see the survey [23] and the references therein). 
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Figure 2: Optimal traffic flow 




If we have previous knowledge about the velocity drift of the distribution of information 
created at the sources (denoted p^) and/or the distribution of information received at the 
destinations (denoted and we assume the Brownian mobility model, then the distribu- 
tion of sources and/or the distribution of the destinations evolves according to the stochastic 
differential equation 

dp+{t) = v+(x, t) dt + a+(x, t) dW^{t) 

and/or dp~(t) = {yi,t) dt + a_{-x,t) dW^ (t). 

where W~^{t) and W~{t) are two independent brownian motion with values in X x Y and 
(T+ := cr+(x, t), cr_ := (T_(x, t) are parameters of the model. 

Assume as in the previous case that we know the initial distribution of the information 
created at the sources. Then by using Ito's lemma, p"*" evolves in time by the Kolmogorov 
Forward Equation 

d d I d"^ 

— p(x, s) = - — [v+(x, s)p(x, s)] + - — [alp{x, s)]. 

for s > 0, with initial condition p(x, 0) = p^(x) 

Equivalently, the initial distribution of the destinations evolves in time by the Kolmogorov 
Forward Equation 

d d _ 1 (9^ 

— p(x, s) = - — [y-{x, s)p(x, s)] + 2^[^-P(^' ^)]- 

for s G [ti,tf), with initial condition p{x,ti) = p^{x). 



5 Optimization on time 

Notice that the minimization problem we solved does not really consider the interaction on 
time because the problem describes the movement of sources and destinations nodes in the 
space and then we solve the static problem at each time. The problem solved at each time 
may not be optimal in the whole period of time considered. 

Another more realistic problem would be to minimize the quantity of nodes used on the 
whole network during a fixed period of time [tj,^/], ^-c, 

/■*/ /■*/ r 

Min / N{D,t)dt= / / r/(x,t)c?x = 

Jti Jti J D 

= f'f |T(x,t)prfx 

Jti J D 

where N{D,t) is the number of active relay nodes in the network D at time t, subject 
to ffT3al) . (fT3b|) . and (TE). 
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For the case we have randomness in the system this problem will be 

Min [ ^ E{N{D,t)}dt= f^f E {|T(x,t)p} rfx 

Jti Jti J D 



subject to fll3ap . (]13b|) . and (TE) since D is compact. From the work of Santambrogio 



page 6) we have the following result: The problem 

Min/ A;(x)|T(x)| such that V ■ T = - z/ 

J D 

is equivalent by duality to the problem of finding 

Min / dk{x,y)d'y such that 7 G n(/i, z/) where 
Jd 

dk{x,y)= inf Lk{uj):= [ k{uj{t))\uj'{t)\dt 

{ui:uj(0) = x. J„ 

In our case /c(x) = |T(x)| then 



Lk{u)= / \T{umWit)\dt 
Jo 

Given that u;(0) = x, and uj{l) = y then by change of variables L^^uj) = JJ|T(x)| dx, and 
as it is independent of u then dk{x,y) = /j^|T(x)| c/x. 

Example.- For the case where we do not have previous knowledge about the velocity drift 
then we just consider the standard Brownian mobility model given by 

dp+{t) = (T+(x,t) dW+{t) 

and / or 

dp-{t) = (T_(x,t) dW'{t). 
where W~^{t) and W~{t) are two independent Brownian motions with values in X x Y. 
Then the previous equations translate into 

d 1 (9^ 

d 19^ 
—p{x, s) = +2^[^' ^)]- 

which have as solution the following equations 
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1 



,2 



e 2t<T 



Now we can replace this solution into Step [2] and Step [3l 

Remark.- Notice that if we suppose that the distribution of the destinations is fixed, as it 
will be the case for aggregation centers of information, then a~ = and then p~{x,t) = p~ 
for all time t. 

6 Conclusions 

We consider a setting in which a spatially distributed set of mobile sources is creating data for 
a spatially distributed set of mobile destinations. We make the assumption that the network 
is massively dense, i.e., there are so many sources, destinations, and nodes, that it is best to 
describe the network in terms of macroscopic parameters, such as their spatial distribution, 
rather than in terms of microscopic parameters, such as their individual placements. We 
focus on a particular physical layer model that is characterized by the following assumptions: 
i) the wireless nodes must only transport the data from the location of the sources to the 
location of the destinations, and do not need to sense the data at the sources, or deliver 
them at the destinations once the data arrive at their physical locations, and ii) the nodes 
have limited bandwidth available to them, but they use it optimally to locally achieve the 
network capacity. 

In this setting, the optimal distribution of nodes induces a traffic flow that resembles the 
electric displacement that will be created if we substitute the sources and destinations with 
positive and negative charges respectively. The analogy between the two settings is very 
tight, and many features of Electromagnetism have a direct interpretation in wireless sensor 
networks. 
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7 Appendix 



Definition.- [Divergence] Let x, y,zhe a system of Cartesian coordinates on a 3-dimensional 
space and let i, j and k be the corresponding basis of unit vectors respectively. The diver- 
gence of a continuously differentiable vector field F = i + Fyj + Fz k is defined to be the 
scalar- valued function: 

_ _ dF^ dFy dF, 

V • F = — - H ^ H -. 

dx dy dz 

An equivalent definition is the following: Given a sequence of surfaces Ak, that all include 
in their interior an arbitrary point {xq, yo, Zq), such that their areas \Ak\ with k, then 

V ■F{xo,yo,zo) = lim -r^ [ F{x,y, z) ■ ndV, 

where n(x) is the unitary normal vector at x. 

We notice that both equivalent definitions can similarly be defined in a 2-dimensional Eu- 
clidean space. 

Definition.- [Curl] Let x, y, z be a system of Cartesian coordinates on a 3-dimensional space 
and let i, j and k be the corresponding basis of unit vectors respectively. The curl of a 
continuosly differentiable vector field F = F^. i -|- j -|- A; is defined to be the vector field 
function: 

^ \ dy dz ) \ dz dx J 
dFy_dFA^ 
dx dy J 

Theorem.- [Divergence Theorem] The divergence theorem states that for any well- 
behaved vector field A defined within a volume V surrounded by the closed surface S the 
relation 

I A-ndS^ I V-AdV 

Js Jv 

holds between the volume integral of the divergence of A and the surface integral of the 
outwardly directed normal component of A. 

The same result holds in a 2-dimensional Euclidean space considering the corresponding 
definition of divergence. 

Theorem.- [Stokes Theorem] The Stokes theorem states that if A is a well-behaved vector 
field, S is an arbitrary open surface and C is the closed curve bounding S, then 



c 



A-d£^ / {VxA)-ndS (19) 
Js 

where n(x) is the unitary normal vector at x. 
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